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The following tables have been compiled to 
facilitate the computation of the orbits of spectro- 
scopic and visual binaries, and small planetary 
objects, such as asteroids and satellites. These 
tables represent a considerable extension of the 
TABLES FOR THE TRUE ANOMALY IN 
ELLIPTIC ORBITS, by Frank Schlesmger and 
Stella Udick, in the Publications of the Allegheny 
Observatory, Volume II, Number 17. Here both 
the true anomaly and the eccentric anomaly are 
given, and the eccentricity ranges from 0.01 to 
1.00. 

The computation involved in the compilation 
of these tables was done on the IBM 709 Data 
Processing System at the Data Processing Center 
of the Agricultural and Mechanical College of 
Texas. The programming technique employed 
was FORTRAN using double precision arithmetic. 

The accuracy of the tables is eight (8) digits. 
The computation was done by the double preci- 
sion method, giving accuracies of ±1 in the 
twelfth digit, and then rounded off was made to 
eight digits. The maximum round off error in 
the tables will be ±1 in the fifth decimal place. 

The line of apsides, PP' in the figure, is the 
major axis of the elliptical path of the minor 
body about the major body. The radius vector, 
FB, is the instantaneous line joining the two 
bodies. Perihelion, P, is the point of closest 
approach of the two bodies. True anomaly, v, 
is the angle at the major body between the line of 
apsides at perihehon and the radius vector. Mean 
anomaly, M, is what this angle, v, would be, were 
the motion of the secondary body uniform. 
Eccentric anomaly, E, is an auxiliary, and may be 



represented as the angle at the center of the ellipse 
between the line of apsides at perihelion and the 
radius of the eccentric circle to the point ob- 
tained by extending to the eccentric circle the per- 
pendicular from the secondary to the major axis. 
The anomalies, E and v, are illustrated in the 
accompanying figure. 

These anomahes are related by the following 
equations : 

E—e sin E=M, (1) 

tan ^ v=J~^ tan jj E, (2) 

2 V 1— e 2 

where e is the eccentricity of the ellipse. Because 
the first equation is transcendental, the task of 
computing E from e and M is one of classic diffi- 
culty and is known as Kepler’s Problem 

Several methods of solving for E were used, 
but all proved unsatisfactory from the standpoint 
of either accuracy or machine time. The first of 
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these methods was a Lagrangian development of 
E in senes (OEUVRES DE LAGRANGE, Vol. 
3, pp. 113-138), 

2 3 

E=M+e sin sin 2M+^ ( 3 * sin 3M 

—3 sin Af)+ • -, (3) 

which, when carried out to the term involving 
e s , gave accuracy only to the second decimal place 
for large e. Several variations of the Newton- 
Raphson Method were tried, and these likewise 
failed to give the desired accuracy. 

The technique finally used was one of utmost 
simplicity, and the actual computation was very 
rapid. (The total time to compute the entries 
in the table and to record these entries in binary- 
coded decimal on magnetic tape was approxi- 
mately 20 hours.) A flow diagram will best 
illustrate our method of solving Kepler’s Equation. 



This procedure was repeated until the absolute 
value of Y was less than 10 -l °. 

To obtain » corresponding to a value of M 
greater than 180°, enter the M - row with a value 
of 360 °—M, and the resulting value of v is then 
subtracted from 360°. For example, for e=0.34, 
M=210°, then 360°— M= 150°, opposite which 
we find fl=360°— 164?08132 = 195?91868. The 
same procedure holds for E. 


Most interpolation will be performed vertically, 
for fractional values of M. The first differences 
are printed in columns adjacent to the values, each 
entry being the difference between the entry cor- 
responding to that value of E or v and the pre- 
ceding. 

The differences between columns are not printed, 
for very seldom does one need to interpolate for 
the third decimal in e. 

Linear interpolation is sufficient m most cases. 
However, the second differences are printed adja- 
cent to the first differences, where agam each entry 
is the difference between the entry corresponding 
to that value of the first difference and the pre- 
ceding To make this correction, the computation 
should proceed as follows : if the fractional part of 
M is less than 0 5, multiply one-half the preceding 
second difference by the fractional part of M, and 
add algebraically to the first difference correction , 
if the fractional part of M is greater than 0 5, 
multiply one-half the corresponding second differ- 
ence by the fractional part of M, and add alge- 
braically to the first difference. 

For most of the table, the error in neglecting the 
second differences will not exceed 10 -2 , except for 
eccentricities between 0.45 and 0.80. However, 
it is suggested that the person using this table 
glance at the second difference column each time, 
to make his own decision. 

The tables were printed directly from magnetic 
tape by means of the IBM 717 Printer, which has 
complete error checking facilities. Hence no typo- 
graphical errors could appear. 

David Larson wrote the program in single pre- 
cision for the IBM 704, but the accuracy was not 
sufficient for our purposes. 

George Taack rewrote and supervised the pro- 
gram in double precision for the IBM 709. 

We wish to express our appreciation to Mr. 
Robert L. Smith, Jr., Head of the Data Processing 
Center of the Texas Engineering Experimental 
Station, for his cooperation and the use of these 
facilities. 
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ECCENTRIC AND TRUE ANOMALY IN ELLIPTIC ORBITS 
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ECCENTRICITY 0.28 0.29 0.30 
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ECCENTRICITY 0.37 0.38 0.39 
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ECCENTRICITY 0.49 0.50 0.51 
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ECCENTRICITY 0.52 0.53 0.54 
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ECCENTRICITY 0.61 0.62 0.63 
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ECCENTRICITY 0.64 0.65 0.66 

M 1ST OIF. 2ND DIF. 1ST DIF. 2ND DIF. 1ST DIF. 2ND DIF 
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CCENTRICITY 0.73 Q.74 0.75 

M 1ST DIF. 2ND DIF. 1ST DIF. 2ND DIF. 1ST DIF. 2ND DIF 
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ECCENTRICITY 0.85 0.86 0.87 

M 1ST DIF. 2ND OIF. 1ST OIF. 2ND DIF. 1ST DIF. 2ND DIF. 
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ECCENTRICITY 0.91 0.92 0.93 
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